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PROCESSOR 


ms  CQMBDUTIOI  OF  TIM  SIRISS  AID  CROeS-SEXn’ICH  DATA 
ZB  niBuiDunirr  flow  aialyszs* 

BiDMthV.  Arrow  \ 

Stanford  mlraralty 

0.  SuiK^y 

Ttw  faliuwlDft  problaa  aroaa  In  tha  couraa  of  a  Largar  ttudy  wnieh  aoxiglrt 
to  axplaln  tba  varlationa  of  lnput*output  ratloa  ovar  tlwa.  (By  tha  In^t-output  » 
ratio  of  Induatry  to  Induatry  ‘J'ls  naaot  tba  ratio  of  that  part  of  tha 
output  of  Industry 'J[  uaad  by  Induatry  'J’to  tba  output  of  Induatry  J,  (aaa  ['ij). 
For  aueb  a  atudy,  tbara  ara  two  typaa  of  data  avallabla.  For  all  yaara,  tbara 
ara  aval labia  (idaally)  outputa  and  final  daaanda  (tba  final  damand  for  an 
Induatry  eonslata  of  all  uaaa  of  Ita  product  otbar  than  in  otbar  Industrlea  or 
Itaalf)  for  all  Industrlaa.  Tba  ^balanea  aguatlona”  of  lnput>output  aaalyala 
(aac  /~3_7  balow)  fora  an  (incoaplata)  ayataa  of  alaultanaoua  ralatlona  wbleb 
any  ba  aatlaatad  by  aona  yaralon  of  tba  nctbod  of  anxlaua  llkallbood  ]for  eoapu- 
tatlonal  raaaona,  tba  a Ingla -aquation  llaltad-lnforaatlon  artbod  la  tba  only 
one  llkaly  to  ba  uaad).'  Bowavar,  for  aoaa  yaara,  wa  bare  additional  Infuraatlon 
In  tba  fora  of  luiovlng  tba  actual  Intarlnduatry 'flows,  which  elaarly  should 
subatantlally  Ineraaaa  tha  accuracy  of  our  aatlantas. 

Tba  alaplast  taehnlqua  la,  of  couraa,  to  aaauat  that  tba  "trua”  input-output' 
ratio  for  any  yaar  for  wbleb  flow  data  ara  avallabla  la  axaetly  aqual  tu  tha 
obaarvad  Input-output  ratio  for  that  yaar.  Tba  asauaptlon  bablnd  this  Is,  bowavar, 
contradictory  to  tha  basic  postulata  that  all  tba  ralatlona  Involvad  ara  valid 
only  up  to  a  stochastic  tarn. 

It  la  tharafura  of  Intarast  to  consldar  aora  axplleltly  tba  Intarlndustry 
flow  aodal  l.apllad  In  tba  usa  of  both  tins  aarlas  and  Intarliulustry  flow  data  for 
Mtlaatli^  tbs  pnraaatars  lnaolwad*r>TtaB  aasiwBi  llkallbood  aatlaataa  for  a  stasia 

•  Tbls  papar  areas  out  of  a  study  eonductad  by  tba  RAID  Corporation.  Z  aa 
Indabtad  to  Ronald  U.  Sbspbard,  now  of  tba  Sandla  Corporation,  for  raising 
tha  problaa. 
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equation  are  then  derived.  In  the  aenae  of  being  expresaed  aa  the  aolutlon 
of  a  aysten  of  sliauitaneoua  equatluna.  Unfortunately,  theae  equatlona  are 
rather  cuabereoae,  and  It  la  hoped  that  aone  reaaonabiy  efficient  oethod  of 

a 

aoeltloB  eaa  be  foead. 


1.  A  Wedel  With  Conatar.t  Input-Output  Coefficianta 

For  expository  reasons,  we  will  start  by  assuxdng  that  the  Input-output 
ratloa  are  cooataat  over  ties.  Mi  eoaaldar  tiM  aatlaatlon  of  a  alofla 
one  of  the  balance  equations  of  the  Input-output  table.  Let 

=  derived  demand  for  coanodlty  0  at  time  l.e.,  total  net  output 
less  final  demand, 

Xji^  ■  output  of  comoodlty  ^  at  time  t  for  J  •  1,  ...,  N, 

Xj^  a  Input  of  commodity  0  to  Industry  at  time 
We  have  the  Identity, 


where  ■  Is  the  number  of  coaaodltles  other  than  coanodlty  0.  The  assu^>tlon  of 
Input^nitput  coefficients  constant  over  time  can  be  written, 

(2) 


where  u^^  la  a  dlaturbanee,  distributed  normally  with  mean  sero.  From  (1)  and 

(2) ,  we  have  the  balance  equation  for  Industry  0, 

(3)  St  •  “‘j  “jt  ♦  "t- 

where 

'  i  "3‘- 

The  disturbance  u^  Is  therefore  normally  distributed  with  naan  0. 
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It  la  aanMarl  th*t  thara  ara  two  typas  of  obaanratlona,  ona  in  which  only 
tba  wariablaa  x.  (J  *0,  M)  ara  obsanrad  and  ona  in  which  tha  varlablaa 

x^^  ara  alao  obaanrad.  Lat  0  be  tha  aat  of  yaara  of  tha  flrat  typM,  V  tha  aat 
of  yaara  of  tha  aaeond  typaj  in  application,  V  containa  only  tha  yaar  1947 

t 

plua  poaaibly  tha  yaar  1939*  Lat 

Vj.t  ■  “jt  “V  ‘ « *• 

x^  ba  tha  coluan  ractor  with  eoaponanta  Xj^  (J  ■  0,  H), 

iL  tba  eoluvi  vaetor  with  coaponants  x  (J  *1,  2M), 

'  Jt 

tba  ooluan  vactor  of  tha  pradateradnad  variablaa, 

IT  tba  ragraaaion  matrix  of  x^  on  x^, 
yr  tba  rafraaaion  matrix  of  x^  on  x^ . 

Tba  raducad  forma  for  can  tbarafora  ba  writtan, 

(5)  x^  ■  77  x^  t 

(6)  \  ^  *t 


whara  tba  ractora  ara  aaeb  diatributad  normally  with  maana  xaro  and  inda- 

pandantly  of  x^. 

Lat  TTj  ba  tha  j'’**  row  of  '77',  tr ^  tha  row  of  tF,  By  aul  atituting  (5) 
and  (6)  into  (1)  and  (2)  and  aquatinf  eoafficianta  of  x^  in  tha  uaual 

(7)  77.*!^  TTj. 

(«)  ®  “j  ^  *  - . “)• 


Prom  tba  dafinition  of 


Vj,t* 


(9) 


Lot  yV.  ba  tba  oovarianoa  matrix  of  ^  that  of  Than  tba  diatribotiona 
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of  and  (for  glvan  x^)  are  given  by 

(10)  f(x^)  ■  k^lAl^  exp  [-(i)(x^  -  x^)*-A.(x^  -  ^  1^^)], 

(11)  g(x^^)  =  ^  exp  [-(4)(x^^  -  W  \)*  ^(\  *^)]» 

where  and  k2  are  constanta.  For  t  in  U,  x^  is  observed;  for  t  in  V,  x^  is 
observed.  Hence,  the  likelihood  function  L  is  obtained  (assusiing  serial 
independence  of  the  disturbances)  by  eultiplylng  together  all  expressions  of 
fora  (10)  with  t  ^  U  with  all  expressions  of  fora  (11)  with  t  7. 

(12)  log  L  =  C  ♦  (T/2)  log  |A|  -  (1)  (x^  -  Tr  »  )*  ./v.  (x.  -  IT  tj 

tfe  U  ^ 

*  (a/2)  log  jAj  -  (i)  (x^  -  ir  1^)*^  (Xg  -  it  x^), 

where  T  and  a  are  nustbers  of  observations  in  U  and  7,  respectively* 

The  ale,  than,  is  to  Baxiaire  log  L  with  respect  to  fT,  ^ 

-/V,  and  A,  subject  to  the  restraints  (7-9)*  The  resulting  values  of 
^2*  ez**  t.he  required  maxi likelihood  estiaates* 


2.  Toe  M^ael  With  Varying  Input-Output  Hatlos 

Ir.  general,  eacn  input-  utput  ratio  Is  presuael  to  vary  linearly  vlth  a 
number  cf  other  va'-'^bles.  Let  be  the  variables  on  which  tue 

Input-output  ratios  o»,  depend.  Any  one  cjefflclent  nay  lepend  on  only  soiae 
of  these  variables;  let  K,  be  the  set  of  variables  on  vhlcno^.  depends.  Por 
slapllclty  of  notation,  let  w  be  tne  constant  1,  and  assuTie  that  K,  contains 


0  for  all  J 


Hg*  5 


(13)  =  E  \- 

(13) and  (2)  can  b«  written, 

(U)  *  (  E  V)  xjt  (J  *  1.  •••.  ")• 

Lat 

*  V  ’jt  <“  =  ® . »'  J  1 . ">• 

For  k  •  0,  thia  dafinltioo  eorraaponda  to  tho  oarlior  one,  sxnca  w  .  ■  1  for 


ot 


all  t.  Alao,  lot 


<»)  ^ . . . «). 

bo  tbo  eoluan  voctor  whoao  coapODonta  aro  (j  a  0,...,  l(E'*'l)), 

x^  bo  tbo  ooltan  roctor  wtaoao  eoapooaots  aro  (j  ■  !,•••,  N(Kr2)). 

Tbo  voctora  x^  and  aro  tbo  oboorvablos  in  tbo  obaonratlono  undor  U  and 
7,  roapoctirolj;  tboj  aro  ao— hat  rodofinod  froa  aoetlon  1.  (14)  oan  bo  written, 

(16)  Xj^  =  E  ^  kj  *  ■k4-J,t  (J  *  1*  •••* 

Ooflna  TT  and  1^,  aa  before,  aa  tho  rofroaalon  aatrleoa  of  x^  and  Z^, 
roapoetlTolj,  on  tbo  rector  of  prodatorainod  rarlablaa  a^.  31oea  aoat  of  the 
coapODonts  of  tho  two  roetora  are  foraad  aa  a  product  of  two  rarlablaa,  there  la 
aa  roaaoo  for  the  rofroaaiono  to  bo  linear}  bowarar,  aa  Andoraon  and  lubin  bare 
■iioaa,  aetiag  aa  if  tbajr  ware  Uoaar  noraal  laada  to  eonalataat  aatlaataa 
(aao  /~3_7*  P*  37^)«  Honeo,  aquatiooa  (5  )  and  (6)  will  atill  bo  regarded 
aa  ralid.  Since  (L)  atlll  bolda,  (7)  la  etill  ralid.  Subatltutlng  (3)  and 
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\ 

(6)  into  (I6)  jrlalida  ttaa  following  annlofiia  of  {6)t 

I 

(17)  -n-j  •  E  «kj  •••» 

Finally »  (15)  iaplles  that, 

(18)  ='^l(k-l)tj  ^  ••*» 

The  likelihood  function  hna  the  eaiae  fora  as  before*  Hence  the 
likelihood  estimates  are  obtained  by  aarlal  ting  (12)  subject  to  the  restraints 

(7) .  (17).  and  (18). 

3*  Restrictions  on  the  CoTarlance  Matrices  of  the  Reduced  Fora  Dlsturbanoes. 

It  can  be  argued  that  the  above  aodel  Iaplles  certain  relations  between 
the  covariance  aatrlces  of  the  reduced  fonts  (5)  and  (6),  l.e.,-A.^  and 
dlnce  (15)  la  an  Ider^tlty. 


MktJ.t  k(k-l)tj,t 


Let  j  be  the  covariance  of  and  ^j*  covariance  of  v^  and  Vy  Then. 

froB  (19) 


‘^IktJ.  Mk'tJ*  •  ®l(k-l)fj, 


.  »(k»-l)tj»  “ 


1.  •••»  K'tli  J»J*  •  1#  *••»  B) 


which  can  also  be  written 

(20)  -  •i-i.J-a  •••»  ■<*♦2);  J  =  ki-1,  ....  M(It2)). 


Slsdlarly.  since  (1)  Is  an  Identity. 

■ 

(21)  »ol  =  ^  (^  =  ■(*♦1))* 

jel  J*!**" 

It  would  be  reasonable  to  Impose  the  conditions  (20)  and  (21)  upon  the 


aaxlBlsatlon  of  the  likelihood  function.  However,  this  would  greatly  coapUcate 
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th«  fon  of  the  estlaation  procedure.  Further,  if  we  wish  to  sssume  that  the 

Interindustry  flows  are  obseirred  only  with  error,  the  relation  (1)  becosMs  a 

stochastic  relation  (if  is  taken  to  refer  to  the  obserred  rather  than  the 

actual  flew),  and  the  arguaent  le^llng  to  (21)  breaks  down.  For  these  reasons, 

* 

we  will  disregard  the  restrictions  on  the  coTariance  Batrlces  of  the  reduced 
fora  disturbances. 

4.  The  Derivation  of  the  Istiaates 
Rewrite  (18)  as 

(22)  (J  -  M^l,  ...,  ■(it2)). 

We  wish  to  aaxiaite  (12)  subject  (7),  (17),  and  (22).  We  use  the  method  of 
Lagrange  aultipliers.  There  are  altogether  M(K^2)-t  1  restrictions,  each  a  vector 
restriction  with  as  aany  coaiponents  as  there  are  predetermined  variables.  To 
each  restriction  we  assign  a  Lagrange  aultiplier,  itself  a  vector.  Let  A 
correspond  to  to  the  J^^  restriction  in  (17),  and  u  ^  to  the 

restrietion  in  (22).  The  Lagrangian  aay  then  be  written, 

(23)  A  s  C-»-(T/2)  log  |^|Ka/2)  log  plj-(l/2)  (x^  -  TT  t^)»yV(Xt  - 

-  (1/2)  Z  ix  ~lr  *,)»_K  (x  -  W  »,) 

se  7 


Let  be  the  aatrix  of  subs  of  cross-products  of  x  and  s,  suaaation  extending 

over  all  observations  in  D, 

M  be  the  aatrix  of  s\a«  of  squares  and  cross-products  of  the  s*s,  svuaoed 
over  observations  in  0, 
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ba  the  matrix  of  sums  of  cross-products  of  x  and  z  over 
be  the  matrix  of  sums  of  squares  and  cross-products  of  the  z*8, 
suBsed  over  V. 

Hov  differentiate  A  successiveljr  with  respect  toTT ,  (k  ■  0, 

J  -  1,  N),  TT'y  and  equate  the  results  to  0. 

(24)  d  k/rr^  =  7\?  (M„  -  7r  f  A*  •  0. 

(25)  dk/d  (11^  -  +  -  0, 

(k  ■  0,  •••«  K{  J  •••»  N)» 

(26)  dh/ ,7?  ^  ^  =0  (J  *  1.  ....  N). 

4 

(27)  dA/^  =Ia.^(m  -  ir  R^j)  -  1/^  »  =  0  (J  =  Ntl.  11(1^2)). 

Here,  each  equation  represents  a  row  vector  of  derivatives.  In  (25),  it  is  under¬ 
stood  that 

(28)  *  0  if  k  ^  ly 
Define  row  vectors, 

(29)  Po  *  -  A*.  1^  t)  (k  e  0,  ...,  K;  J  =  1 . K), 


(30)  Hj  =  /^  ♦  •  (  j  =  1,  M).  Rj  =  ♦  (J  =  8(1^2)), 


and  matrices. 


(31) 

>o 

Pi 

[*°  1 
“i 

P  = 

• 

/■(li-li 

p  R  • 

• 

• 

• 

h{if2) 

Then  (24)  and  (25)  may  be  written. 


(32)  .  n  «„)  .  p. 


(26)  arid  (27)  b«C( 
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(33)  X(i  -  tri  )  =  R. 

Jx  •• 

fro.  (29)  and  (30). 

(34)  kj  *  *Sl(k»l)»J  “  •‘kj  • 

Dafina 

(35)  Qoq  •  ■  “'kj*  %k4j,I(ktl)f-J  *  *  0  ^o**  9  ^ 

R(k4>l)4>J  (q  a  1,  ...,  Kl^R);  k  <  0.  ....  I;  J  »  1.  ....  H). 


Lat  Q  ba  tha  aatrlx  with  alaaanta  (p  ■  0.  ....  ■(l-tl)i  q  ■  1.  ....  M(K^2)). 
Lat  Q  ba  tha  oolavi  raetor  with  ^  «  1.  *  *°^kj  (^  *  ^*  ***•  ^>  J  ~  •••* 


M). 


Than  (34)  can  ba  wrlttan. 

(36)  P  s  Q  B  -t-  Q  . 

o 


If  wa  ooabina  (32)  and  (36). 


(37)  -A-(II^  -  TTM,,)  =  Q  R  ♦  Q  . 

Iota  that  Q  and  Q  dapand  onijr  on  tha  In  affaet.  tha  unknown  Lagranfa 

paraaatars  ^  (^  *  ***»  ^  (J  *  •••»  M(l^2))  hava  baan 

tranaferaad  into  tha  unknown  aatrix  R  and  vector  P^.  (33)  and  (37)  are  tha 

raaulta  of  tha  differentiation  of  the  Lagrangiao  (23)  with  respect  to  the 
ragrasslon  aatricas  tr.  TT.  Row  diffarantiata  with  respect  to 

(3a)  -  o  (k  =  o.  ....  ij  j  =  i.  ....  a). 

idith  the  aid  of  (29)  and  (30).  (36)  aay  ba  rewritten 

(^’^  ^RkfJ  ("j”*"  **o^’  ■  (•*  ■  **  J  *  •••• 
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We  can  elioinate  Pg  as  follows:  In  (37) »  oultiply  on  the  left  by  Q* 
where  ■  _A.”^  •  Let 

(LO)  y  -  (Q*  Z  Q)"^  • 

Note  that  y  is  a  scalar.  Then, 

(a)  Pg  =  y  (Q*  -  Q*.  ‘  0*  T,  Q  R)- 

s 

finally,  we  diffsreritiate  A  with  respect  to  the  aatrices  .A.  , 
the  inwerses  of  the  cowariance  matrices  of  reduced  form  disturbances  in  the 
equations  in  the  observations  of  U  and  V,  respectively.  Let 
as  in  Anderson  and  Subin  (see  n.  2), 

(W)  r  =■  Ixx  -  "-X.  tr’  -  ^  ic„  »  Tm,.  Tr'  i 

(43)  r  =  W*  -  TT  i  .  fr  R  . 

'We  have  then  to  solve  equations  (33 )»  (37),  (39)  i  (42),  and  (43)  for  the 
unknown s_A ,  R»  If  »  7T,  and  TT,  where  P^  is  eliminated  with  the  aid  of 

(40)  and  (41).  Actually,  we  are  only  interested  in  Q  (note  that  the  matrix  ^ 
is  completely  defined  by  the  vector  Q)  which  involves  the  structural  coefficients 
to  be  estimated.  These  equations  seem  considerably  more  difficult  to  handle  thar 

e 

the  corresponding  ones  of  Anderson  and  Eubin. 
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